Astronomy Sz Astrophysics manuscript no. ferrigno'lE1145 


© ESO 2008 


February 2, 2008 





INTEGRAL observation of the accreting pulsar 1E1145. 1-6141 

Carlo Ferrigno^, Alberto Segreto^, Teresa Mineo^, Andrea Santangelo^, and Riidiger Staubert^ 



lASF-Pa, via Ugo la Malfa 153, 90146 Palermo Italy 

lAAT, Abt. Astronomie, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 



o 
o 

(N 

Q 

oo 



O 



O 



Received February 2, 2008; accepted 



Abstract 



Aims. We analyze 1050 ks of INTEGRAL data of the high mass X-ray binary pulsar IE 1145.1-6141 to study its 
properties over a long time baseline, from June 2003 to June 2004, with wide spectral coverage. 

Methods. We study three high luminosity episodes, two of them at the system apoastron, three brightening with lower 
intensity, two at the periastron, and one extended period of intermediate luminosity spanning one orbital cycle. We 
perform timing analysis to determine the pulse period and pulse profiles at different energy ranges. We also analyze 
the broad band phase average spectrum of different luminosity states and perform phase resolved spectroscopy for the 
first flare. 

Results. From the timing analysis, we find a pulse period of ~ 297 s around MJD 53 000 with a significant scatter 
around the mean value. From the spectral analysis we find that the source emission can be described by an absorbed 
bremsstrahlung model in which the electron temperature varies between ~25 and '~37keV, without any correlation to 
luminosity, and the intrinsic absorbing column is constantly of the order of 10^"^ cm~^. Phase resolved spectral analysis 
evidences a different temperature of the plasma in the ascending and descending edges of the pulse during the first fiare. 
This justifies the pulse maximum shift by ~ 0.4 phase units between 20 and 100 keV observed in the pulse profiles. 
Conclusions. The comparison with the previous period measurements reveals that the source is currently spinning- down, 
in contrast to the long term secular trend observed so far indicating that at least a temporary accretion disk is formed. 
The study of the spectral property variations with respect to time and spin phase suggests the presence of two emitting 
components at different temperatures whose relative intensity varies with time. 
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1. Introduction 

X-ray Binary Pulsars (XRBPs) were discovered more than 
, 30 ye ars ago with the pio neer observation of Cen X- 
■ 3 by iGiacconi et al.l (|l97lh . Although the basic mech- 
I anisms of the pulsed emission were quickly und erstood 
! (jPringle &: ReesI 119721 : [Davidson fc Ostrikeil [l973l) . many 
' puzzling aspects still persist. In particular, their wide band 
spectral behavior has not been completely explained on the 
'• basis of a self consistent physical model; spectra are suc- 
cessfully described by phenomenological models such as a 
power law modified at high-energy by an exponential cut- 
off and at low energy by intrinsic photo-el ectric absorp- 
tion plus a thermal compo nent at ~ 0.1 keV (jCoburn et al.l 
I2002t ldi Salvo et al.ll200i ). XRBPs are powered by the ac- 
cretion of ionized gas from a nearby companion, an O or B 
star (M > SMq) in High Mass X-ray Binaries (HMXBs), a 
later than type A star or a white dwarf {M < Mq) in Low 
Mass X-ray Binaries (LMXBs) onto the surface of a black 
hole or neutron star (NS). If the compact object is highly 
magnetized, with surface field as high as 10^^ G, the plasma 
is threaded at the magnetospheric boundary, and then 
funneled onto the m agnetic poles fPringl e fc ReesI Il972t 
[Davidson fc Ostrikeil 1 1973.) along the field lines; pulsation 
is generated because of the nonalignment between the mag- 
netic and the rotational axes. In HMXBs the plasma can be 
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captured directly from the slow and dense equatorial stellar 
wind of a Be star, or from the radiation driven wind of an 
early type super-giant. In brighter systems, accretion can 
be also powered by Roche-lobe overflow via an accretion 
disk. 

The systems with Be companions (often with apprecia- 
ble orbital eccentricity) are normally quiescent but some- 
times they flare when the NS passes through the compan- 
ion's disk. During these episodes the luminosity may in- 
crease by a large factor, a temporary accretion disk may 
form, and the NS generally spins up. The bright HMXBs 
accreting from a stable accretion disk via Roche lobe over- 
flow have rather short spin and orbital periods, nearly cir- 
cular orbits, and can exhibit torque reversal episodes su- 
perim posed on a spin-u p trend as observed for Cen X- 
3 (Bildsten et all 119971) . Most wind-fed HMXBs with a 
super-giant companion, have long pulse periods (hundreds 
of seconds), moderate and rather stable luminosity {Lx — 
2Q35-37 gj-g g-i-j^ and can spin-up or spin-down depending 
on the transfer of angular momentum from the wind to the 
magnetosphere. The pulse period evolution is then charac- 
terized by red-noise fluctuations due to th e irregular ac- 
cretio n flow as in the case of Vela X-1 (jBildsten et al.l 
Il997f ). Simulations have shown that a prograde or retro- 
grade accretion disk, which is not necessarily a true vis- 
cous disk, may form and disrupt on time scales of 10^ s, 
causing changes in t he sign of the pulse period derivative 
([Blondin et al.|[T"990r ). as observed in the controversial case 
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of the system with unknown companion OAO 1657-415 
(lBavka]|ll997l :ll nam temporary accretion 

disk is also suggested to explain the short episodes of spin- 
up, accompanied by luminosity enhancements, in the wind- 
fed binary GX 301-2 (iKoh et alJ Il997t llnam fc Bavkal 
[2001 . 

Thanks to the Galactic plane monitoring, performed 
with the INTEGRAL satellite, a new class of HMXBs 
containing a neutron star and super-giant (SG) donors 
has been identified: the Super-giant Fast X-ray Transients 
(SF XT), characterized by short outbursts of X-ray emis- 
sion (|Negueruela et al.ll2006l:rSguera et al.ll2005ll2006[ ). The 
transient behavior of these sources is explained by the mo- 
tion il l wide orbits ( 10 stellar radii) through a clumpy SG 
wind (|in't Zandll2005 iNegueruela et al.ll2006t) . or to the re- 
current pass age of the NS thro ugh a denser equatorial wind 
component (jSidoH et al.ll2007D . 

In this paper, we report on the INTEGRAL obser- 
vations of one of the less known HMXBs located very 
close (170 to the exten sively studied 4U 1145-61 (see 
e.g. IStevens et a"Llll997D . IE 1145.1-6141. Optical spec- 
troscopy of its B2 lae s uper-giant companion yi elds a dis- 
tance of (8.5 ± 1.5) kpc (jDensham fc Charle£| [l982). X-ray 
spectral analysis was first performed using the data in 
the energy range 0.5-3 keV from th e imaging propor tional 
counter of the Einstein observatory ()Lamb et al.lll980D find- 
ing a column density of A^h = (3 ± 2) x lO^^cm"^. This 
value was confirmed, using the proportional counter ar- 
ray of_thfi_^RXTE_satelliteJ^ the energy range 2-60 keV, 
by iRav fc Chakrabartvl (|2002l ) who modeled the spectrum 
using an absorbed cut-off power law with index F — 1.4, 
cut-off energy Ec = 6.4 keV and e-folding energy Ef = 
18 keV. They also measured the orbital parameters of the 
binary system finding an eccentric (e = 0.20) orbit with 
period P = 14.365 days, asini = 99.4 It-sec, periastron 
epoch To = 51008.1 MJD, and longitude of periastron 
u! = —52°. They proposed, on the basis of previous ob- 
servations, a se cular spin-up trend of ~ 1 0^^ s/s between 
1975 and 2000 (|Rav fc Chakrabartvl [200l) . The NS emis- 
sion was characterized by a pulsed fraction of ~50% in 
the 4-20 keV energy band and an X-ray luminosity of 
-Wff_erg/s (Hutchings et al.lll98ll; ICrebenev et allll992t 
iRav fc Chakrabartv 200 j) . Such a low luminosity is in- 
consistent with Roche lobe overflow for this class of bi- 
naries and indicates that the NS is almost certainly ac- 
creting from the companion wind. The pulse shape has 
been studied firstly using the ART-B telescope on board of 
the GRANAT observa tory in the 3-60 keV energy range by 
iGrebenev et al.l (|1992D who found evidence of a notch be- 
tween the maximum and the minimum of the pulse; this fea- 
ture appeared also in the p ulse shape obtained with RXTE 
(|Rav fc Chakrabartvl[200a ). 

Our paper is aimed to study the spectral and 
timing properties of IE 1145.1-6141 between March 
2003 and June 2004 using the data from European 
Space Agency's International Gamma-Ray Astronomy 
Laboratory {INTEGRAL) and is organized as follows: in 
Sect. [5] we describe the observations, in Sect. [3] we illus- 
trate the scientific results which are discussed in Sect. |4l 
and in Sect. [5] we draw our conclusions. 



2. Observation and data analysis 

INTEGRAL, launched in October 2002, carries 3 co-aligned 
coded mask telescopes. 

- IBIS (Imager on B oard the INTEGRAL SateUite; 
lUbertini et all I2003D . which allows for 12 arcmm 
FWHM angular resolution imaging in the energy range 
from 15 keV to 600 keV with energy resolution of 8% at 
100 keV and 8° fully coded field of view. It is composed 
by a low energy Cd Te detector (15-600 keV), ISGRI 
(jLebrun et a ll I2003D. and by a Csl layer (175keV- 
lOMeV), PICsIT (iLabanti et al.l[200l . 

- SPJ_ (spectrometer on INTEGRAL; IVedrenne et~all 
120031 ) ■ sensitive from 20keV to 8MeV with an angu- 
lar resolution of 2.5° and an energy resolution of a 0.2% 
at 1.3 MeV; the fully coded field of view is 16°. 

- JEM -X (Joint European X-ray Monitor; iLund et al.l 
|2003| ) that includes two independent units JEM-Xl and 
JEM-X2, sensitive from 3 keV to 34 keV, with an angu- 
lar resolution of 1 arcmin and an energy resolution of 
13% at lOkeV; its fully coded field of view is 4.8°. 

INTEGRAL observations are divided in short periods 
(of ^ 2000s each), called science windows (SCWs), during 
which the telescope maintains the same pointing. From one 
science window to another, the pointing follows a dithering 
or a survey strategy as chosen in the observing program 
(iWinkler et al.l[200l . 

To reduce and analyze JEM-X data we used the Off-line 
Science Analysis (OSA) software version 5.1; for ISGRI data 
we ex ploited also the software d e scribe d in iFerrigno et al.l 
()2007D and ISegreto fc Ferrignol (|2007f) . For the spectral 
analysis we used the standard XSPEC package version 12.2. 

In this paper we used the publicly available INTEGRAL 
data in which IE 1145.1-6141 was visible before June 2004. 
Thanks to the ISGRI large field of view it has been possible 
to observe the source for a total of 1.08 Ms between June 
2003 and June 2004, with a nearly continuous coverage of 
2 binary orbits in 2004, and a sparse coverage of two and a 
half orbits in 2003. We limited our analysis to ISGRI and 
JEM-X data because the source is too weak to be studied 
with SPL As summarized in Table HI JEM-X detected the 
source only in a sub-set of the available SCWs, because of 
its smaller field of view with respect to ISGRI. 

As it can be seen from the light curves of Fig. [1] we 
detected three flares of about the same intensity (labeled 
HI, H2 and H3, where "H" is an abbreviation for High 
luminosity), the first two occurred in correspondence to 
the NS apoastron, the last while the NS was moving away 
from the companion. Apoastron flares of this source have 
been observed also by th e Burst Alert Telesc ope on board 
of Swift and reported bv lCorbet et all (|2007D . 

At MJD 53 144 (day 1600 in the scale of Fig. [T]) the 
source was so weak that it was pointless to perform any 
analysis. In the following orbital cycle, it exhibited weak 
and irregular luminosity increases, none of which was in- 
tense enough to allow a dedicated analysis. We then treated 
these data as a unique observation named EI, meaning that 
the source had an Extended Intermediate luminosity level. 
We divided the rest of the observations into three intervals 
labeled LI, L2 and L3 (where "L" is an abbreviation for 
Low luminosity), during which the source exhibited less 
pronounced brightening. In the interval selection we ex- 
cluded the neighboring SCWs with few source counts to 
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improve the S/N and we separated L3 from EI since it has 
the appearance of a low luminosity flare occurring near the 
system periastron, such as LI, rather than an extension of 
EI. We finally note that the L2 episode occurred while the 
NS was approaching the companion. 

Th e H3 episode has been reported by iBodaghee et al] 
(|2004[ ). who stated that a flare of this intensity has never 
been mentioned in the literature for IE 1145.1-6141. But it 
should be noted that the source did present brightening of 
the same intensity during the episodes called here HI and 
H2. 

To perform the timing analysis, we transformed the 
event arrival times for each time interval listed in Table [T] 
to the Solar System Barycenter and then to the Binary 
System Barycenter using t he orbital parameters found by 
iRav fc Chakrabartvl (|2002f ). We determined the pulse pe- 
riod with a version of the phase-s hift method described ex- 
tensively in iFerrigno et al.l (|2007f ) and briefly summarized 
here. Using the ISGRI data, which guarantee the best sig- 
nal to noise ratio (S/N) and time coverage, we find a ten- 
tative pulse period from the Fourier analysis of the light 
curve. Then, using this value, we accumulate a pulse pro- 
file in each science window and determine the phase of the 
pulse maximum. Finally, from a parabolic fit of the pulse 
maximum phase as a function of time, we find the optimal 
pulse period and an upper limit on its first derivative. 

Unless differently stated, the errors in the paper are 
quoted at la confidence level. 

3. Results 

3.1. Timing analysis 

We measure the pulse period in the time intervals reported 
in Table [2l while during LI, L2 and L3 the statistics is 
not sufficient for an independent determination. The pulse 
period derivative is always too small to be directly mea- 
sured. In particular, during EI we find the 2a upper limit 

P < 4 X lO^^s/s. Thus, during this orbital cycle the pe- 
riod variations do not exceed the long term trend. 

We plot in Fig. [5] our results to gether with the pe- 
riod a nd the spin-up trend reported bv lRav fc Chakrabartvl 
(|2002f) . We note that the INTEGRAL values are signif- 
icantly higher than the extrapolation of the s ecular evo- 
lution measured by iRav fc Chakrabartvl (|2002D . To verify 
the occurrence of a torque reversal, we computed the aver- 
age value of our period determinations (296.695 ± 0.002 s) 
and of the last four RXTE results (296.572 ± 0.001s): 
the difference of 0.123 ± 0.003 s in five years evidences an 
average spin-down of ~ 7 x 10^^°s/s in contrast to the 
previous spin-up trend of —1 x 10~^s/s. The source has 
then begun to reduce its average angular speed on a time 
scale of years even though fluctuations from the trend are 
present. In fact, at MJD 52824.0 (H2) we found a period 
of (296.62 ± 0.02) s, compatible with the RXTE determina- 
tions, but at more than 3a fro m the ot her values we mea- 
sured. If the trend proposed bv lRav fc~ Chakrabartv (2002) 
is correct a torque reversal would have taken place between 
RXTE and INTEGRAL observations around MJD 52 000; 
but we cannot exclude that the source was spinning down 
before MJD 45 000, spinning up until MJD 50 500 and then 
spinning up again. Unfortunately, the poor determination 
of the pulse period by EINSTEIN at MJD 44000 does not 
allow to verify this second hypothesis. 
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Figure 1. Light curve of IE 1145-6141 during the 
INTEGRAL observations. In each panel wc plot, from top 
to bottom, the data of JEM-X (4-25 keV) and ISGRI (22- 
71keV), (from http://isd.c.unige.ch/?Data+sources) 
and the distance between the stars of the binary system 
expressed in l ight seconds as computed fr om the ephemeris 
measured by iRav fc Chakrabartvl (2002) . The time inter- 
vals in which we performed the analysis are indicated in 
the figure between vertical solid lines and labeled with the 
corresponding name. 



Using the values of Table [21 we produced the pulse pro- 
files of the ISGRI data in the energy bands 15-35, 35- 
60, and 60-llOkeV setting arbitrarily the zero phase at 
the pulse minimum in the lowest energy range. The pulse 
profiles (Fig. [3]) show a similar shape over all the observ- 
ing periods. They are dominated by a single broad peak 
which moves to earlier phase for increasing energy, and 
present also a hint for a second peak at phase ~ 0.7, 
which is reminiscent of the notch between the pulse maxi- 



mum and minimum report ed bv iGrebenev et al.l (jl992) and 
IRav fc Chakrabartvl (|2002D . However, some differences arise 
at higher energy where the main peak has an increasing 
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Table 1. Summary of the INTEGRAL observations of IE 1145.1-6141. 



Name 


c;j_„„j_ lA/T TT\| 

btart [IVIJJJJ 


r^na [MJJJJ 


Exposure [ksj 


Mean Count Rate [cts/s] 


Ho-oU key ) 








ISGRI 


JEM-X 


ISGRI 


JEM-X 


[erg s ] 


Ll 


52 788.8947 


52 791.0241 


64.1 


10.9 


5.9±0.1 


3.2±0.1 


4.1x10'"'' 


HI 


52 795.6284 


52 796.4538 


29.2 


13.7 


14.4±0.2 


6.2±0.1 


1.1x10^^ 


L2 


52 814.7897 


52 815.5495 


49.6 


15.5 


5.8±0.2 


4.3±0.5 


2.5x10^^ 


H2 


52 823.7611 


52 824.9043 


39.0 




14.7±0.2 






EI 


53 146.8561 


53 162.3828 


703.5 




3.53±0.05 






L3 


53 162.6755 


53 163.3679 


42.1 




7.9±0.2 






H3 


53 165.3538 


53 169.1147 


153.4 


2.2 


12.2±0.1 


8.35±0.03 


1.0x10^*^ 



assuming a distance of 8.5 kpc. 



Table 2. Pulse periods (P) and the upper limits on its first derivatives (P). 



name reference time* [MJD] 


P[s] 




P 


[s/s] 


HI 52 796.00566 


296.78±0.04 


< 7 X 10"" 


H2 52 824.25522 


296.62±0.02 


< 6.5 X lO"'' 


EI 53 154.80607 


296.698±0.003 


< 4 X 10"^ 


H3 53 167.38765 


296.69±0.01 


< 1. 


8 X 10"^ 



*At the reference time the pulse profile has a minimum in the 15-35 keV energy range. 
**Upper hmits are given at the 2a confidence leveL 
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Figure 2. Long-term evolution of the pulse period of 
IE 1145.1-6141. The diamonds after MJD 52 000 represent 
the results of our analysis. The crosses are the previous de- 
terminations, and the solid line is the spin-up trend derived 
bv lRav fc Chakrabartvl (|2002l) . 



sharpness, if the time intervals are put in the sequence EI, 
H3, HI, and H2. 

To investigate the evolution of the pulse profile with 
energy, we focus on the HI time interval when JEM- 
X exposure is long enough to accumulate pulse profiles 
with an acceptable S/N below 20keV. The pulse profiles 
are modeled with a Gaussian plus a constant (reduced 
between 0.8 and 1.4 for 35 degrees of freedom) as shown 
in Fig. m We find that the width of the Gaussian remains 
constant with an average value of 0.21 ± 0.05 while its cen- 
troid shifts at earlier phase for increasing energy: between 
2 and 8keV the pulse profile peaks at pulse phase 0.5 ±0.2, 
between 20 and SOkeV at pulse phase 0.4 ± 0.1, and above 
SOkeV at pulse phase 0.2 ± 0.2. The linear fit of the phase 




Pulse Phase 



Figure 3. Background subtracted pulse profiles of 
IE 1145.1-6141 for different periods in three energy bands. 



of the centroid position x as function of energy E gives: 

TP 

X = (0.55 ± 0.10) - (4.4 ± 2.4) x 10"=^- (1) 

[keVJ 

Repeating the same procedure for the intervals H2, EI and 
H3 for the ISGRI data alone, due to the bad quality or lack 
of JEM-X data, we find that the relation of equation ([T]) 
is always verified. However, the limited energy range pre- 
vents an independent determination of the centroid evolu- 
tion with energy. 

3.2. Spectral analysis 

In each observing period we perform broad band phase 
averaged spectral analysis using ISGRI data and, when- 
ever available, JEM-X data. We investigate several spec- 
tral models and eventually find satisfactory fits for cut- 
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Figure 4. Background subtracted pulse profiles of 
IE 1145.1-6141 during HI. Below 20keV the data are from 
JEM-X , above from ISGRI . In each panel, the solid line 
is the best fit pulse profile obtained adding a Gaussian to 
a constant; the vertical dashed line is the position of the 
Gaussian centroid. 



off power law {E"'* exp[—E / Ef]) plus photoelectric ab- 
sorption, i mplemented as PHABS in XSPEC using the cross 
sections by iBalucinska-Church fc McCammonl (1993), and 
bremsstrahlung also attenuated by absorption. But a dis- 
crimination between these models on a statistical base is 
not possible because all fits give acceptable values (see 
Table O. 

The cut-off power law model gives best fit values with 
large uncertainties in most of the observing periods prob- 
ably due to a strong correlation between the spectral in- 
dex and the other parameters. It is then difficult to find 
any significant trend of the parameters with luminosity or 
time; nonetheless some conclusions can be drawn: the low- 
est value of iVn (7.5±2.5 x 10^^ cm^^) is significantly higher 



than both the Galactic column (TVh = 1.4 x 10 cm 
iDickev fc Lockmaiil Il990f ) and the intrinsic c olumn mea- 
sured by EINSTEIN ((3 ± 2) x lO^^cm'^, iLamb et al.l 
|198C1() : the measured spectral index o f the power law (1- 
1.5) is typical of several HMXBs (Co burn et al.ll2002[ ). 

In Fig. [5l we show the spectra in the low and high 
luminosity states LI and HI together with the high en- 



Ray & Chakrabarty (2002) 




Energy channel [keV] 



Figure 5. Broad band unfolded energy spectra for the high 
state HI and the low state LI. Below ~20 keV data are from 
JEM-X, above this value from ISGRI. The model in both 
cases is absorb ed bremsstrahlung. The da shed line is the 
model used by iRav fc Chakrabartvl ()2002f ) with arbitrary 
normalization. 



ergy cut-off model used by iRav fc Chakrabartvl (|2002f ) 
(E-'y exp[-{E - Ec)/Ef] for E > E^ and E-"/ for E < 
Ec). We note that in our data the source emission is consid- 
erably harder, as confirmed by the higher value of the fold- 
ing energy Ef — 28^3 keV when the s pectra during HI and 
H3 are described by the same model of lRav fc Chakrabartvl 
()2002D . Fitting the spectra of the other observing periods 
with this function, we note unreasonably large uncertainties 
on the parameters together with a not better than the 
values of Table [31 It is then evident that the relatively low 
S /N of the INTEGRAL observations of this source prevents 
the use of this phenomenological model, which has a redun- 
dant number of free parameters. We finally note that, since 
the source brightness in our observations is not greater than 
in the past, we can exclude that this evolution is luminosity 
related. 

From the best fit values of the bremsstrahlung model 
plotted in Fig. [6l we note that the plasma temperature 
is variable even if not correlated with the source luminos- 
ity. Instead, iVu is always higher than the Galactic column 
and does not present any significant trend. We also ver- 
ify that fixing the absorbing column at its average value, 
the determinations of the plasma temperature do not vary 
significantly. 

3.3. Phase resolved spectral analysis 

During the time interval HI the statistics allows to perform 
broad band phase resolved spectroscopy combining data 
from both ISGRI and JEM-X. In the other time intervals 
the source is either too weak or too off-axis for JEM-X to 
collect enough counts, and then to significantly constrain 
the continuum parameters. 

In this analysis we adopt the absorbed bremsstrahlung 
model, since it gives a good description of the source emis- 
sion with the lowest number of parameters. We divide the 
pulse in five equally spaced phase intervals and fit the spec- 
tra obtaining a reduced ^ ^ for ^100 d.o.f. and the the 
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Figure 6. Best fit parameters of the absorbed 
bremsstrahlung model as function of the source lumi- 
nosity. In the left panel we plot the plasma temperature 
and in the right panel the Nn of the intrinsic photoelectric 
absorber. Errors are at 90% c.l. for one parameter of 
interest. 



best fit parameters plotted in Fig. [T] The value of A^h is 
approximately constant lO^^cm"^), while the plasma 
temperature is higher at pulse phase 0-0.4 (ascending edge 
of the pulse) than at pulse phase 0.4-1 (descending edge). 
Averaging the fitted values, we find fcro-o.4 = (40 ± 2) keV 
and fcro.4-1 = (29 ± 1) keV corresponding to a variation 
of the temperature with a significance of Scr. Fixing the 
value of iVn at 10^^ cm~^ (the average over the spin pe- 
riod) we obtain temperature determinations compatible to 
the previous case, and a slight improvement in the signif- 
icance of the variation (S.Str). The higher temperature in 
the rising part of the main peak explains naturally its shift 
with energy since at increasing energy the hotter emission 
is prominent. 

The phase resolved spectral analysis in the other time 
intervals does not give a statistically significant variation of 
the plasma temperature. 



1 1 1 1 1 1 1 
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Figure 7. Phase resolved spectral analysis results. Upper 
panel: the pulse profile during HI in the 20-60 keV energy 
range with the adopted phase intervals. Lower panel: the 
best fit parameters of the adopted model (photo-electric 
absorbed bremsstrahlung) as a function of the pulse phase 
during HI. From top to bottom we plot the hydrogen col- 
umn density A^h, the plasma temperature kT, when iVu is 
a free parameter and when A^h is held fixed at lO^'^cm"^. 



4. Discussion 

We obtain a good description of the source spectral energy 
distribution using the absorbed bremsstrahlung model; this 
does not mean, alone, that this is the real emission mech- 
anism. To make a consistency check, we compute from the 
expression o f the thermal br emsstrahlung emissivity, equa- 
tion (87) of iLongaiil (|1999D . the quantity nlV, where Ue 
is the electron density and V is the volume of the emitting 
region for a pure hydrogen plasma emitting at the observed 
temperature with a power compatible to the luminosity of 
the source. Assuming for the accretion column a radius of 
100 m and estimating the mass accretion rate from the X- 
ray luminosity, it is possible to derive an order of magni- 
tude estimation for the electron density which implies a 
Compton optical depth of order unity both in the vertical 
and in the horizontal directions of the accretion column. 
Therefore, it is not possible to ensure that Compton scat- 
tering is negligible, but, on the other hand, it is not clear 



how much it distorts the underlying spectral energy distri- 

bu tion. 

iBecker fc Wolfj (|2007D recently showed that the ob- 
served spectrum of XRBPs is well described by thermal and 
bulk Comptonization of the optically thin bremsstrahlung 
and cyclotron emission from the accretion column. Using 
the measured luminosity and distance of the source, we 
tested their model on the spectrum with best S/N, i.e. 
during HI. The fit was very good evidencing that all 
the emission comes from re-processing of an underlying 
bremsstrahlung with an electron temperature of ~ 13keV 
without any need for intrinsic absorption. Unfortunately, 
it was not possible to constrain all the model parameters, 
especially the magnetic field, due to the absence of absorp- 
tion features. However, we note that the spectral energy 
distribution is very similar to absorbed bremsstrahlung and 
then, in data with relatively low S/N, the two models are 
statistically equivalent. 
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Table 3. Summary of the results of spectral fits with two models: cut-off power law and thermal bremsstrahlung, both 
modified by photoelectric absorption. The fluxes are computed directly from the source counts and are thus independent 
from the adopted absorption; for comparison the flux of the Crab in the 20-50 keV energy band is lO^^ergs"^ cm~^. The 
errors are at 90% confidence level for one parameter of interest (Xmin + 2.7). 



parameter 


Ll 


HI 


L2 


H2 


EI 


L3 


H3 


Cut-off power law 
















Nh [10^2 
Photon index 


20 ±8 
1.5 ± 0.4 


7.0 ±2.5 

1.1 ± 0.1 


1.5 ±0.5 


N/A 
1.4 ±0.2 


N/A 
1.8 ±0.3 


N/A 
1.1 ±0.8 


12 ±7 
1.2 ±0.2 


Ef [keV] 
xVd.o.f. 


218/197 


25.5 ±2.5 
232/220 


39tf3 
240/205 


28^1 
36/43 


or + 18 

64/51 


25tf 
74/65 


197/195 


Bremsstrahlung 
Nh [10^2 cm-2] 
kT [keV] 


15 ±4 
28 ±2 


10.4 ± 1.5 
33.3 ± 1.5 


15 ±9 
37 ±4 


N/A 
28.5 ± 1.0 


N/A 
25.1 ± 1.3 


N/A 
35 ±4 


13 ±5 
32 ± 1 


xVd.o.f. 


220/198 


240/221 


240/206 


36/44 


68/52 


62/44 


198/196 


flux(5-15keV)* 
flux(20-50keV)* 


2.9 ±0.2 

2.9 ± 0.1 
_i 


7.47 ±0.15 
8.2 ± 0.1 


2.3 ±0.5 
2.9 ± 0.1 


N/A 
8.8 ±0.1 


N/A 
1.96 ±0.04 


N/A 
4.0 ±0.2 


6.6 ±0.3 
7.29 ±0.07 



* in unity of 10 ^"^ergs ^ cm 



The values we find for the plasma temperature (higher 
than for the Becker and Wolff model) and for the ab- 
sorbing column depth should then be regarded as a useful 
parametrization of the spectral energy distribution, rather 
than a measurement of the plasma characteristics. The 
strong absorption could be due to an emission turn-off at 
low energy, but could also be partially explained by the 
presence of optically thick gas in the stellar wind e.g. a bow 
shock formed by the X-ray pulsar radiation (Blondi n et al.l 
Il990l) . Unfortunately, only observations with a better spec- 
tral resolutions below ~ 10 keV would detect the absorption 
edges and then constrain significantly the column value. 

More insight into the structure of the emission region 
can be achieved from the pulse phase resolved spectroscopy. 
The constant N-^ at any phase indicates that the low en- 
ergy part of the spectrum is governed by a process isotropi- 
cally distributed around the NS or located far from it, while 
the bimodal temperature distribution and the presence of 
a notch in some pulse profiles (Fig. [3]) suggest the pres- 
ence of two components at different temperatures in the 
magnetosphere of IE 1145.1-6141. If their relative inten- 
sity is variable with time, the different values of the plasma 
temperature in the phase averaged spectra between 25 and 
40keV and the time modulation of the pulse shape would 
be naturally justified. Indeed, all the analyzed spectra can 
be fitted with a linear combination of two bremsstrahlung 
model with fixed temperatures of 19 and 40keV. 

During three flares, the source reached 80-90 mCrab in 
the 20-50 keV energy band. During three other periods it 
reached 30-40 mCrab, but for one binary orbit (EI) it pre- 
sented only irregular and weak luminosity variation with 
an average flux of 20mCrab. Two of the most intense flares 
occurred at the apoastron, the other at about two thirds of 
the maximum distance. At the periastron the NS reached 
at most one third of the maximum luminosity. 

This behavior rules out the hypothesis of spherically 
symmetric wind, and of a combination of disk and wind 
proposed to explai n the periodic fl ares in IGR J11215-595 2 
(|Sidoli et al.ll200l or GX 301-2 (iPravdo fc G hoshl [200lh . 
A complex wind environment with clumps or filaments of 
higher density plasma is a more natural explanation of its 
irregular luminosity increases, but our measurements do not 



allow a good spectral constraint of the absorbing column, 
thus this hypothesis cannot be tested spectroscopically. 

A clumpy stellar wind is invoked also to explai n the light 
curves of the super-giarit fast X-ray transients (jin't Zandl 
l2005t iLevder et al.l l2007l ) characterized by short irregular 
flares lasting from a few ks to a fraction of a day. The 
NS in these systems is thought to orbit at ^ 10 stellar radii 
while in persistent systems the separa tion is supposed to be 
smaller () Walter fc Zurita Herasl[2007[) . We can estimate for 
IE 1145.1-6141, on the basis of the spect ral type, a mass of 
^ 10 Mq and a stellar radius of ~ 5 Rq (|Habets fc Heintz^ 
|1981D . From the measured orbital parameters, we derive 
that the NS lies at about 10 stellar radii from the compan- 
ion, with a distance modulation of 20% due to the orbit 
eccentricity. 

Although it is a persistent system, it has some charac- 
teristics similar to the SFXT: in particular the component 
separation and the flaring activity. However, compared to 
these objects, the flares are longer (from a significant frac- 
tion of a day to about four days) and have a smaller rela- 
tive (4-5) intensity with respect to the persistent flux level, 
which, in turn, is two orders of magnitude larger. Supposing 
that the mechanism to capture the plasma has similar effi- 
ciency in all systems, we infer that the wind is on average 
denser, while the clumps are less compact than in SFXT, 
making the system an intermediate case between transient 
and persistent wind-fed HMXBs with super-giant compan- 
ion. 



Unfortunately, in our case, an estimate of the mass spec- 
trum of the clumps from the flare duration is not straight- 
forward since the torque history is characterized by local 
variations, and then the formation of a temporary accretion 
disk is probable, which prevents the plasma from directly 
falling onto the NS surface. Wind clumpyness appears a 
common feature in binary early type stellar wind suggest- 
ing further investigations to clarify whether this can be due 
to the interaction with the NS or it is an intrinsic charac- 
teristic of the wind. 
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5. Conclusions 

We confirm that IE 1145.1-6141 is a persistent low lumi- 
nosity binary system but we find interesting peculiarities 
not previously noticed. 

The source was always detectable, when in the instru- 
ment field of view, and showed a complex light curve which 
strongly suggests a clumpy stellar wind. 

The timing analysis shows that the pulse period of the 
source is not constant (at 98% confidence level). These vari- 
ations of the torque are expected in wind- fed systems, in ad- 
dition we find that the source reversed its long term torque. 
In fact, the averaged of the measured pulse periods is higher 
than the average of the RXTE observati o n (J une 1997- 
January 2000) while iRav fc Chakrabartvl ()2002l) found a 
previous spin-up trend. Unfortunately, the lack of continu- 
ous monitoring prevents to study in more detail the torque 
history and then the accretion model. 

The absorbed bremsstrahlung model allows an incisive 
parametrization of the source spectral characteristics, even 
though it cannot be regarded as a coherent physical model, 
since the Compton optical depth of the accretion column 
is of order of unity. We find that the plasma temperature 
varies with no correlation to the source luminosity while 
the absorbing column remains approximately constant. A 
comparison with the model adopted in t he previous RXTE 
observations (|Rav fc Chakrabartvll2003 ) indicates a signifi- 
cant hardening of the spectrum which cannot be explained 
by a different luminosity level. 

Phase resolved spectroscopy shows that the absorbing 
column depth is approximately constant, and the plasma 
temperature in the first half of the pulse profile is (11 ± 
2) keV higher than in the second half. This justifies the 
pulse maximum shift by ^ 0.4 phase units between 20 and 
lOOkeV. To describe the emission of IE 1145.1-6141, we 
then propose the presence of two spectral components at 
different temperatures and, possibly, a layer of absorbing 
material surrounding isotropically the NS or located far 
from it. The variation of the components' relative intensity 
with time would also explain the oscillations of the phase 
averaged temperature that we observe. 
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